o 

(N 

0^ 






c^ 



o 



rn 

O 

o 



X 



Polarization dependence of coherent phonon generation and detection in 
highly-aligned single-walled carbon nanotubes 

L. G. Booshehri,i'2 c. L. Pint,2'3.4 q jy Sanders,^ L. Ren,i'2 C. Sun/- 2 E. H. Haroz,i'2 
J.-H. Kini,6 K.-J. Yee,6 Y.-S. Lim/ R. H. Hauge,^'^ C. J. Stanton,^ and J. Konoi-^'^.g 

'Department of Electrical and Computer Engineering, Rice University, Houston, Texas 77005, USA 

^The Richard E. Smalley Institute for Nanoscale Science and Technology, Rice University, Houston, Texas 77005 

'^Department of Physics and Astronomy, Rice University, Houston, Texas 77005, USA 

^Department of Chemistry, Rice University, Houston, Texas 77005, USA 

^Department of Physics, University of Florida, Box 118440, Gainesville, Florida 32611-8440 

^Department of Physics, Chungnam National University, Daejeon, 305-764, Republic of Korea 

''Department of Applied Physics, Konkuk University, Chungju, Chungbuk, 380-701, Republic of Korea 

(Dated: July 20, 2010) 

We have investigated the polarization dependence of the generation and detection of radial breath- 
ing mode (RBM) coherent phonons (CP) in highly- aligned single-walled carbon nanotubes. Using 
polarization-dependent pump-probe differential-transmission spectroscopy, we measured RBM CPs 
as a function of angle for two different geometries. In Type I geometry, the pump and probe po- 
larizations were fixed, and the sample orientation was rotated, whereas, in Type II geometry, the 
probe polarization and sample orientation were fixed, and the pump polarization was rotated. In 
both geometries, we observed a very nearly complete quenching of the RBM CPs when the pump 
polarization was perpendicular to the nanotubes. For both Type I and II geometries, we have de- 
veloped a microscopic theoretical model to simulate CP generation and detection as a function of 
polarization angle and found that the CP signal decreases as the angle goes from 0° (parallel to 
the tube) to 90° (perpendicular to the tube). We compare theory with experiment in detail for 
RBM CPs created by pumping at the £^44 optical transition in an ensemble of single- walled carbon 
nanotubes with a diameter distribution centered around 3 nm, taking into account realistic band 
structure and imperfect nanotube alignment in the sample. 

PACS numbers: 78.67.Ch, 63.22.-(-m, 73.22.-f, 78.67.-n 



I. INTRODUCTION 



The one-dimensionality of single-walled carbon nan- 
otubes (SWNTs) is attractive from both fundamental 
and applied points of view, where the ID confinenient 
of electrons and phonons results in unique anisotropic 
electric, magnetic, mechanical, and optical propertiesi^ 
Individualized SWNTs, both single-tube and in ensem- 
ble samples, have shown anisotropy with polarized Ra- 
man scattering and optical absorption measurements 
where maximum signals result when the nanotube axis is 
aligned parallel to the polarization of incident light<^"— 
Additionally, due to strong anisotropic magnetic sus- 
ceptibilities, both semiconducting and metallic SWNTs 
align well within an external magnetic field, and with the 
added properties of the Aharonov-Bohm effect, the elec- 
tronic band structure of SWNTs respond anisotropically 
with the strength of a tube-threading magnetic fluxi^""— 
Such optical and magnetic anisotropy is also expected 
with bulk samples, but detailed measurements showing 
extreme anisotropy have been lackingJi 

Here we investigate anisotropic optical and vibrational 
properties of a bulk film of highly-aligned SWNTs us- 
ing polarization-dependent coherent phonon (CP) spec- 
troscopy. CP spectroscopy is an ultrafast pump-probe 
technique that complements CW Raman spectroscopy, 
and although both techniques provide information about 
electron-phonon coupling, CP spectroscopy avoids the 



common disadvantages of Raman spectroscopy that in- 
clude detection of Rayleigh scattering and photolumi- 
nescence and the broadening and blending of peak 
features.^- Such advantages are useful when investigat- 
ing SWNTs, where a large majority of samples include 
ensembles dispersed in various environments and their 
optical properties are obscured by the collection of vary- 
ing species of nanotubes. 

Recent CP studies on SWNTs have produced direct 
observation of CP oscillations of both the radial breath- 
ing mode (RBM) and G-band phonons, in addition to 
their phase information and dephasing times j^^Ti^I Fur- 
thermore, when pulse-shaping techniques are combined 
with CP spectroscopy, predesigned trains of femtosec- 
ond pulses selectively excite RBM CPs of a specific 
chirality, avoiding inhomogeneous broadening from the 
ensemble32. However, it is important to note that pre- 
vious CP studies investigated randomly-aligned SWNT 
samples, and as the quasi- ID nature of SWNTs leads to 
optical anisotropy that is dominant in the polarization 
dependence of PL, absorption, and Raman scattering, 
CP studies on aligned SWNTs are necessary. 

Below, we present results of a detailed experimental 
and theoretical study of CPs in highly- aligned SWNT 
thin films and found a very strong polarization anisotropy 
of the RBM as a function of angle. In particular, we ob- 
served a very nearly complete quenching of the RBM 
when the optical polarization is perpendicular to the 
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FIG. 1: (color online) Absorption spectra of aligned single- 
walled carbon nanotubes in the (a) near-infrared to visible 
and (b) far-infrared (or terahertz) ranges for parallel (Ay) 
and perpendicular (Ax) polarization. 



tubes. We have developed a theoretical model to un- 
derstand this extreme anisotropy, including band struc- 
ture, interband optical transition elements, ultrafast car- 
rier and phonon dynamics, and imperfect nanotube align- 
ment in the sample. Fitting our results to theory, we also 
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calculated the nematic order parameter, 5' 
degree of alignment of SWNTs in the sample. 



i.e., the 



II. SAMPLES AND EXPERIMENTAL 
METHODS 

We investigated CPs through degenerate pump-probe 
spectroscopy measurements on highly-aligned SWNT 
thin films transferred onto sapphire substrates. Pat- 
terned, vertically aligned SWNT arrays grown by chemi- 
cal vapor deposition were subsequently etched with H2O 
vapor to allow transffer to sapphire via a dry contact 
transfer printing technique, which forms horizontally- 
aligned SWNT thin films.— »SS, The resulting film pro- 
duces aligned nanotubes of the same length, with a di- 
ameter distribution centered around 3 nmi^ Figure [TJa) 
shows polarization-dependent absorption spectra of a 
typical aligned sample in the visible to near-infrared 
range, while Fig. [Ijb) shows absorption spectra in the 
far-infrared or terahertz (THz) range. The polarization 
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FIG. 2: (color online) Scanning electron microscopy image of 
aligned single-walled carbon nanotubes, and the two exper- 
imental configurations employed in the current pump-probe 
spectroscopy work, (a) Type I: pump and probe polariza- 
tions are fixed and sample orientation is rotated, (b) Type 
II: probe and sample orientations are fixed and pump polar- 
ization is rotated. 



anisotropy is obvious in both (a) and (b), but extreme 
in the THz regime (b), as there is virtually zero ab- 
sorption when the sample is perpendicular to the THz 
polarizationJi Therefore, with such an aligned sample, 
our polarization measurements can be extended to in- 
clude the sample as an additional rotation parameter, 
compared to our earlier study of CP polarization depen- 
dence in randomly-aligned SWNTsi^ 

Using a mode-locked Ti:Sapphire laser with ^80 fs 
pulse width that is shorter than the period of excited 
CP oscillations and ~40 mW average pump power, the 
laser was tuned to central wavelength of 850 nm (1.46 eV) 
to predominantly excite the £'44 interband transitions of 
the SWNTs in the sample. A shaker-delay system and 
balance detector was employed for fast-scan, real-time 
observation of CPs. Extraction of CP oscillations from 
raw pump-probe time-domain signal were performed as 
previously described.-^ As depicted in Fig. [21 two types 
of polarization measurements were investigated. Type I 
configuration maintained the same polarization for the 
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FIG. 3: (color online) Experimental difTerential transmission 
data in the time domain showing coherent phonon oscillations 
of the radial breathing mode for different polarization angles 
in (a) Type I and (b) Type II configurations (see Fig. [T]). The 
traces are vertically offset for clarity. 



pump and probe, while the alignment axis of the sample 
was rotated. Type II configuration maintained the same 
orientation for the probe and sample, while the pump 
polarization was rotated. For Type II measurements, a 
half- wave plate provided 90-degree rotation of the pump. 
All measurements were performed at room temperature. 



III. EXPERIMENTAL RESULTS 

Results of our polarization-dependent ultrafast pump- 
probe differential transmission measurements on the 
aligned SWNT film for both Type I and Type II ori- 
entations are shown in Fig. [3] Here, it is seen that the 
typical differential transmission amplitude of the REM 
CP oscillations is on the order of 10~^, with a CP de- 
cay time of ~1.5 ps. This short decay time, compared 
to our earlier studies on individually-suspended SWNTs 
in solutionji2ii2ii2^ is expected with our sample of bun- 
dled SWNTs. As the polarization angle is rotated, we 
see a strong polarization anisotropy of the REM CPs as 
a function of angle, where the strongest oscillations are 
observed at 0° while the signal is very nearly completely 
quenched at 90° for both Type I and Type II geometries. 

Figure |4] shows CP spectra for different polarization 
angles for Type II configuration. To produce these spec- 
tra in the frequency domain, we calculated the Fourier 
transform (FT) of the time-domain CP oscillations. The 
frequencies of the CP oscillations have a wide range, from 
50 to 200 cm~^, which is consistent with the large di- 
ameter distribution of nanotubes within the sample.^^ 
The polarization anisotropy is clearly observed in the CP 
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FIG. 4; (color online) Coherent phonon spectra for different 
polarization angles in Type II configuration obtained through 
Fourier transform of the time-domain data in Fig. ^h). The 
traces are vertically offset for clarity. 
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FIG. 5: (color online) Spectrally integrated coherent phonon 
intensity as a function of angle, measured in (a) Type I and 
(b) Type II configurations. 



spectra. Figures [5ja) andlSjb) plot the integrated CP in- 
tensity of the FT for all excited nanotubes as a function 
oi 9 [in Fig. [SJa)] and cj) [in Fig.ISJb)], where 9 {(f) is the 
angle of rotation for the Type I (Type II) configuration. 



IV. THEORY 

We have developed a microscopic theory for the gener- 
ation of coherent phonons in single-walled carbon nan- 
otubes and their detection in coherent phonon spec- 
troscopy experiments. Our microscopic theory is de- 
scribed in detail in Ref. |22| . so we only summarize the 
main points here and indicate how our earlier work has 
been extended to include polarization-dependent coher- 
ent phonon spectroscopy. 

We treat the tt and tt* electronic states in {n,m) car- 
bon nanotubes using a third-nearest-neighbor extended 
tight binding (ETB) formalism developed by Porezag et 
al^ for carbon compounds. Using a density functional 
based parametrization, Porezag et al. derived analytic 
expressions for the Hamiltonian and overlap matrix ele- 
ments that depend only on the C-C bond lengths. Using 
the ETB formalism, we obtain tight-binding wave func- 
tions and electronic energy levels Eg^ (k) where s = c,v 
labels the conduction and valence bands, /i labels the cut- 
ting lines, and k is the one-dimensional nanotube Bril- 
louin zoncii By exploiting the screw symmetry of the 
nanotube, we can block-diagonalize the electronic Hamil- 
tonian and overlap matrices into 2x2 blocks, one for each 
cutting line, as described in Appendix A of Ref. 123 . 

We treat nanotube lattice dynamics using a seven- 
parameter valence-force-field model described in Ap- 
pendix B of Ref. 22. Exploiting the nanotube screw sym- 
metry, the dynamical matrix can be block-diagonalized 
into 6x6 blocks, one for each cutting line, which can 
be solved for the phonon displacement vectors and dis- 
persion relations uj'i^{q). In the the dispersion relations, 
/3 = I,--- ,6 labels the phonon modes, v is an angu- 
lar momentum quantum number, and q is the phonon 
wave vector in the one-dimensional nanotube Brillouin 
zone. Following the work of Jiang et al."^^ and Lobo and 
MartinS)2S we include four types of force field potentials, 
i.e., the bond stretching, in-plane bond bending, out-of- 
plane bond bending, and bond twisting potentials. Our 
force-field potential energies are invariant under rigid ro- 
tations and translations (force constant sum rule), and, 
as a result, our phonon model correctly predicts the dis- 
persion relation for the long-wavelength fiexure modes-^^ 

The electron-phonon interaction is treated in a second- 
quantized formalism where the electron-phonon inter- 
action matrix elements are evaluated in Appendix C 
of Ref. [2^. In calculating electron-phonon matrix ele- 
ments, we use 2pj, graphene atomic wave functions and 
screened atomic potentials obtained from an ah initio 
calculation.^- 

We obtain equations of motion for coherent phonon 
amplitudes using the Heisenberg equations as described 
by Kuznetsov and Stanton,'^*' We assume that the opti- 
cal pulse and photoexcited carrier distributions are dis- 
tributed uniformly over the nanotube. In this case, only 
the V = q = Q phonon modes are excited. To excite 
RBM coherent phonons, the laser pulse must be short 
in comparison with the RBM phonon oscillation period. 



For coherent RBM phonons, the coherent phonon ampli- 
tude is proportional to the tube diameter D{t). The tube 
diameter, being proportional to the coherent phonon am- 
plitude, satisfies a driven oscillator equation 



d''D{t) 



+ Lo^D{t) ^S{t) 



(1) 



where lo is the angular frequency of the v = q = Q RBM 
phonon mode. The driving function S{t) for RBM co- 
herent phonons is given by 



S{t) 



sfj,k 



Sfl 



{k)[U{k,t)~f^^{k)] (2) 



where Ms^{k) is the driving function kernal for RBM co- 



herent phonons defined in Eq. (12) of Ref.l22land fs/iik, t) 
and fgn_{k) are the time-dependent and initial equilibrium 
carrier distribution functions, respectively. 

From Eq. ([2]) we see that the driving function depends 
on the photoexcited carrier distribution functions. We 
treat photoexcitation of carriers in a Boltzmann equa- 
tion formalism and obtain the photogeneration rate us- 
ing Fermi's golden rule. In Ref. i22 we only considered 
linearly-polarized laser pulses with the electric polariza- 
tion vector parallel to the nanotube axis. For parallel po- 
larization, optical transitions only occur between states 
with the same cutting line index /i, and the resulting 
photogeneration rate is given by Eq. (13) in Ref. 123. 

In the present work, we consider linearly-polarized 
pump and probe beams in which the electric polariza- 
tion vector makes a finite angle with the nanotube axis. 
If e is the unit electric polarization vector for the pump, 
then the photogeneration rate is now given by 



dfs^dk) 



dt 
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TT^e^ u{t) ff^ 
[fs'^.'{k,t) - fs^.{k,t)) s(^AE^^^:{k) - hi^ 



s' fl' 



(3) 



where AE^^, (k) = \Esf,{k) - Es'f,'ik)\ are the k- 
dependent transition energies, fku is the pump photon 
energy, u{t) is the time-dependent energy density of the 
pump pulse, e is the electronic charge, ttiq is the free 
electron mass, and ng is the index of refraction in the 
surrounding medium. To account for spectral broaden- 
ing of the laser pulses, the delta function in Eq. ([3|) was 
replaced with a Lorentzian lineshape^ 



S{AE - huj) 



Tp/i27r) 



[AE - h^Y + (rp/2) 



(4) 



In Eq. (j3]), P]^^ (k) are the dipole- allowed optical ma- 
trix elements between the initial and final states where fi 
and ii' are initial and final state cutting lines. For polar- 
ization parallel to the tube axis (taken to be z), ^' — fi 
and z ■ Pgg^ik) = Pgg/{k) with the right hand side being 
defined in Eq. (14) of Ref. [22|. For linear polarization 



parallel to the x axis, /^' = /i ± 1 and in the notation of 
Ref. u2, we have 



x-P^f^\k)- 



lY.c:,{s',ti±i,k) 



rJ 

(5) 
Similarly, for linear polarization along the y axis, we have 



^,p,.,±i^f^^^^^^^,p,,,±i^,^^-^ 



':fl,fldzl / 



(6) 



In Eqs. ([5]) and (jS]), the atomic dipole matrix element 
vectors (which can be evaluated analytically) are given 
by 

M(r', rJ) = y dr ^l,^{v - Rr-o) V V5rj(r - Rrj) (7) 

where the 2pz orbitals (fr3 are defined in Eq. (C3) of 
Ref. m. 

In our experiments, a probe pulse is used to mea- 
sure the time- varying absorption coefficient. In our the- 
ory, the time- varying optical properties measured by the 
probe pulse are obtained from the imaginary part of the 
dielectric function. For a linearly polarized probe pulse, 
we have 



£2{hw,t) 



8^^e 



2„2 



SS' flfl' 






(8) 



where At = 7r((it/2)^ is the cross-sectional area of 
the tube and df is the equilibrium nanotube diameter. 
The dirac delta function is replaced by a broadened 
Lorentzian of the form shown in Eq. ^. The photon 
energy of the probe pulse photons is hw, and e is the 
probe unit polarization vector. In our experiments, the 
photon energy fvjj is the same for both pump and probe, 
while the polarization vector e may be different for pump 
and probe. To obtain the computed coherent phonon 
spectrum, we take the power spectrum of the computed 
differential transmission signal after background subtrac- 
tion using the Lomb periodogram algorithm described in 
Ref. [32|. 



V. COMPARISON OF THEORY AND 

EXPERIMENT 

To compare experiment with theory, we performed 
simulations of polarization-dependent CP spectroscopy 
on a (38,0) zigzag nanotube. This is a mod-2 semicon- 
ducting tube with a diameter of 3.01 nm. Our sample 
contains an ensemble of nanotubes with diameters cen- 
tered around 3 nm,^ and we expect that (38,0) tubes 
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FIG. 6: (color online) (a) Theoretical bandstructure for (38,0) 
nanotubes showing the strong £44 and £77 transitions for a 
polarization vector parallel to the tube. The wavevector k is 
expressed in units of tt/T where T = 4.31 A is the length of 
the translational unit cell, (b) Computed absorption spectra 
for polarization angles varying from 0° (parallel to tube) to 
90° (perpendicular to tube), (c) Coherent phonon spectra as 
a function of pump energy and pump polarization angle (j) for 
RBM coherent phonons {ui — 10.91 meV). Coherent phonon 
spectra are Fourier transforms of computed time dependent 
differential transmission for Type II pump-probe experiments, 
(d) Integrated power (black dots) obtained by taking the area 
under the computed £44 peaks in lower left panel. We fit (red 
curve) the theoretical calculations to Acos^{4') where A and 
p = 4.017 are fitting parameters. 



will contribute strongly to the CP signal for the ensem- 
ble since (i) the (38,0) tube has a diameter in the center of 
the measured diameter distribution and (ii) mod-2 zigzag 
tubes tend to have very strong intrinsic CP signals J^^i^ 

Our theoretical results for the (38,0) nanotube are 
shown in Fig. [Bj Our computed electronic tt band struc- 
ture for the (38,0) tube is shown in Fig.IHl^a). The bands 
are doubly degenerate with two distinct cutting lines cor- 
responding to each band. The tt valence bands have neg- 
ative energy, and the tt* conduction bands have posi- 
tive energy. For z-polarized light, the allowed -E44 and 
Ejr transitions (selection rule A/z = 0) giving rise to the 



strongest CP signals are indicated by vertical arrows. 

The absorption coefficient for (38,0) nanotubes as a 
function of linearly-polarized photon energy is shown in 
Fig.m^b) for polarization angles varying from 0° (parallel 
to tube) to 90° (perpendicular to tube). 

For the (38,0) nanotube, the RBM phonons have a 
computed energy hui = 10.91 meV and an oscillation 
period of 379 fs. We simulate the generation of RBM 
coherent phonons by pumping with a 50 fs pump pulse, 
which is much shorter than the RBM oscillation period. 
Polarization-dependent CP spectra for RBM coherent 
phonons are shown in Fig. [BKc) for pump polarization 
angles (/) ranging from 0° to 90°. In these simulations, 
the probe polarization is kept fixed parallel to the tube 
axis (Type II geometry). As can be seen in the figure, the 
CP signal is maximum when the pump is polarized paral- 
lel to the tube. As the pump photon energy is varied, we 
excite RBM coherent phonons by successively photoex- 
citing carriers in different bands. The strongest RBM 
CP signals are obtained by pumping at the E44 and E^j 
transitions. As the pump polarization angle increases, 
the RBM CP signal decreases until it is finally quenched 
when the pump polarization is perpendicular to the tube. 

By taking the area under the RBM CP signal curves in 
the vicinity of the computed -E44 transition energy near 
1.48 eV, we obtain the integrated power as a function of 
pump polarization angle shown as black dots in Fig.[5Jd). 
The integrated CP power can be well fit with a fitting 
function of the form AcosP{9), where A and p = 4.017 
are fitting parameters. 

The polarization dependence of the RBM CP inte- 
grated power in Fig. ^d) can be understood by examin- 
ing Fig. [71 In Fig.[7fb) we plot the coherent RBM diame- 
ter oscillations as a function of the pump polarization an- 
gle 9. Fitting the theoretical results to AcosP{6), we ob- 
tain p = 2.11 as indicated in the figure. For the RBM di- 
ameter oscillations we get p ~ 2 which makes sense. The 
driving function S{t) is proportional to the photogener- 
ated carrier density as seen in Eq. ^ and this in turn is 
proportional to the squared optical matrix element which 
has a cos^(6') angular dependence'^. In Fig. [TJ^a) we plot 
the amplitude of the time-dependent differential trans- 
mission oscillations measured by the probe beam after 
background subtraction. For the differential transmission 
signal, we get p ~ 2 for Type II and p ~ 4 for Type I. 
This also makes sense. In Type II experiments the probe 
polarization is parallel to the tube axis and the pump po- 
larization angle 6 is varied. In this case the signal should 
be proportional to the amplitude of the diameter oscil- 
lations. In Type I experiments on the other hand, the 
pump and probe polarizations are parallel and we should 
get an extra factor of cos^(0) to account for anisotropy 
of the probe beam absorption. In our CP spectroscopy 
measurements, we extract the power spectrum of the dif- 
ferential transmission oscillations. The CP power spec- 
trum should be proportional to the square of the am- 
plitude of the differential transmission signal which we 
anticipate will give us a cos^(0) polarization dependence 
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FIG. 7: (color online) In (b) we plot the amplitude of RBM 
coherent phonon oscillations as a function of the pump polar- 
ization angle (6* or 0) for a (38,0) nanotube photoexcited by 
a 50 fs pump at the theoretical E44, transition. In (a) we plot 
the amplitude of the resulting differential transmission signal 
for Type I and II experiments as a function of angles 6 and (j> 
respectively. Our results are fit to Acos^{6) [A cos^ {(!>)] where 
numerical values of the best fit p are indicated in the figure. 



for Type II CP spectroscopy experiments and a cos^(6') 
dependence for Type I CP spectroscopy experiments. 

The experimental integrated CP power is obtained by 
taking the Fourier transform power spectrum of the time- 
dependent differential transmission data shown in Fig. [3] 
and then computing the area under the £^44 peak in the 
resulting spectrum. The results for our Type I and Type 
II experiments are shown in Fig. |S1 where the experi- 
mental data points are plotted as downward-pointing red 
triangles and the corresponding theoretical predictions 
are shown as upward-pointing black triangles. In over- 
laying the experimental and theoretical data points, we 
subtracted a background from the experimental data to 
obtain complete quenching at 90°. We then rescaled ex- 
perimental and theoretical data (both in arbitrary units) 
so that the integrated power at 0° is equal to unity. We 
then fit experimental and theoretical data to functions 
of the form AcosP{e) [AcosP((/))]. The best fit functions 
for our experimental and theoretical results are shown in 
Fig. [5] as solid red and dashed black lines, respectively. 

For the Type II experiments, where the probe po- 
larization is fixed parallel to the average tube orienta- 
tion, wc get decent agreement between theory and ex- 
periment with best fit parameters p = 4.017 for theory 
and p = 3.630 for experiment. For the Type I exper- 
iments, where pump and probe polarization vectors are 
parallel to each other, there is a discrepancy between the- 
ory and experiment. As can be seen in Fig.|SlJa), the best 
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FIG. 8: (color online) Integrated RBM coherent phonon 
power as a function of (a) 6 for Type I and (b) (j> for Type 
II experiments pumping close to the E44, transition. Experi- 
mental data points are red downward pointing triangles fit by 
a solid red line. Theoretical data points for a (38,0) nanotube 
are black upward pointing triangles and are fit by a black 
dashed line. The fitting functions are of the form Acos^{9) 
and A cos*'((/)) where the numerical values of the best fit p are 
indicated in the figure. The experimental curves have under- 
gone background subtraction to obtain complete quenching 
at an angle of 90° and then rescaled so the integrated power 
at an angle of 0° is equal to unity. 



fit parameters are p = 8.308 for theory and p = 4.435 for 
experiment. The experimental fits in Fig. |8] would imply 
that the CP intensity scales roughly as cos^(0) indepen- 
dently of the probe polarization. This seems unlikely 
given the anisotropy of the absorption coefficient in car- 
bon nanotubes. In fact, from our theory, we would expect 
p ~ 4 for Type II and p ^ 8 for Type I. 

We believe that the cause of this discrepancy is most 
likely due to misalignment effects. Our sample consists 
of an ensemble of nanotubes lying horizontally on a sap- 
phire substrate. While the tubes are highly-aligned, they 
are not perfectly aligned. To consider the effects of mis- 
alignment on the experimentally-measured integrated CP 
intensity in a Type I experiment, we assume that the tube 
alignment angles on the transferred film are described by 
a Gaussian distribution function with a small standard 
deviation A9. If the angle between the pump polariza- 
tion vector and the ensemble averaged tube axis is 0, then 
the angles •& between the pump polarization vector and 
the axes of each tube in the ensemble are described by a 
Gaussian distribution 



P{^,0,A9) 



1 



27r(A6l) 



exp 



2(A6i)2 



(9) 



If the integrated RBM coherent phonon power for each 
nanotube in the ensemble is given by Acos^(^), the en- 
semble averaged integrated power Icp{0,A6) is obtained 



FIG. 9: (color online) Effects of tube misalignment on inte- 
grated coherent phonon power for an ensemble of nanotubes 
with tube axis orientation angles following a Gaussian distri- 
bution, (a) Type I fitting function Acos^{6) with p = 8 and 
(b) Type II fitting function Acqs'' {(j)) with p — 4. 



by taking the ensemble average over -d 

/oo 
dd P{d,9,A0) cosP(t9) (10) 
-00 



In Eq. (|T0|) . we extended the integration limits on i? to 
infinity in the limit of small A9. If p is a positive in- 
teger, Icp{9,A6) can be found analytically. Otherwise, 
it can be evaluated numerically. In the case of Type II 
experiments, the ensemble averaged integrated power is 
obtained by replacing 6 with <j> in Eqs. (JH) and (fTO|). 

The effects of nanotube misalignment on the integrated 
CP power fitting function for an ensemble of Gaussian 
misaligned tubes is illustrated in Fig. [9l assuming p = 
8 for Type I and p — 4 for Type II experiments. As 
the standard deviations Ad [A(p] increase, quenching of 
the CP intensity is reduced and Icp{0,A6) [Icp{<j), Acfi)] 
become flatter. 

Using this model, we are able to get reasonable fits to 
the experimentally measured CP intensity in both Type I 
and Type II experiments. We assume a fitting function of 
theiormAcosP{e + A0)+B [AcosP{(j) + A(j)) + B], wheie 
A and B are background subtraction and rescaling pa- 
rameters and A6 [At/)] is a random Gaussian distributed 
misalignment angle (standard deviation A9 [Ac/)]). In 
our fitting procedure, we set p = 8 for Type I and p — 4 
for Type II and use the same standard deviation for the 
tube misalignment angles in fitting both Type I and II 
data. The results of our fitting procedure are shown 
in Fig. [TUl where the best fits for Type I and Type 
II geometries are shown as dashed black and solid red 
lines, respectively. Our best fit standard deviation is 
A9 — A0 — 18.7°, which implies that the tube align- 
ment angles on the sapphire substrate are 0° ± 9.35°. 
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FIG. 10: (color online) Experimental integrated RBM coher- 
ent phonon power for the E44 transition as a function of 8 for 
Type I (black upward pointing triangles) and for Type II 
(red downward pointing triangles) experiments. Type I exper- 
iments are fit to A cos^{6 + A6) + B where p = 8 (black dashed 
line) and Type II experiments are fit to A cos^{(f>+A(t>) + B for 
p = 4 (solid red line). A and B are background subtraction 
and rescaling parameters and the standard deviations for the 
random tube axis misalignment A6 and Acp axe restricted to 
be the same for both Type 1 and II. 



For small A0 (measured in radians) the nematic order 
parameter is 5" = cxp(-2 {AO)^) = 0.81. 

Although a value of 5* = 0.81 indicates a strongly 
aligned sample, it is not perfectly aligned {S = 1). This 
is unlike previous work in the THz regime with these 
highly- aligned samples, where the nematic order param- 
eter was calculated to be exactly S* = 1.^^ It is clear that 
there is a wavelength dependence of the nematic order 
parameter, but we believe this can be explained quali- 
tatively: any slight misalignment in the sample would 
go undetected in the THz regime, as the wavelength is 
much larger than our visible wavelengths that can de- 
tect the misalignments. This could explain the differ- 
ence in calculated nematic order parameter between our 



CP measurements and previous THz measurements, but 
nonetheless, regardless of the calculated differences with 
wavelength, our calculated nematic order parameter is 
still quite large and we can confidently confirm the high 
degree of alignment of the sample with our CP measure- 
ments. 



VI. CONCLUSION 

In summary, we investigated the polarization 
anisotropy of coherent phonon dynamics in highly- 
aligned single-walled carbon nanotubes and measured 
RBM coherent phonons as a function of polarization 
angle. We saw a very nearly complete quenching of 
the RBM for both geometries and extended our results 
to determine the degree of alignment of the sample. 
Comparing our results with theory, we performed 
simulations of polarization-dependent CP spectroscopy 
on a (38,0) zigzag nanotubes and also found a decrease 
in CP signal as optical polarization varies from parallel 
to perpendicular to the nanotube axis. Using those 
simulated results, we also theoretically determined a 
cos^(6') dependence for Type 1 CP spectroscopy experi- 
ments and a cos*((/)) polarization dependence for Type 11 
CP spectroscopy experiments. Including misalignment 
effects to our fitting, we finally determined the nematic 
order parameter of our sample to be 0.81. 
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